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Abstract 

If three right-handed neutrinos are added to the Standard Model, then, for the 
three known generations, there are six quarks and six leptons. It is then natural 
to assume that the symmetry considerations that have been applied to the quark 
matrices are also valid for the lepton mass matrices. Under this assumption, the 
solar and atmospheric neutrino data can be used to determine the individual neutrino 
masses. Three minima have been found, using the x 2 fit) and, from these minima, 
it is determined that the mass of the lightest neutrino is 1.3 x 10 -3 eV, that of the 
next heavier neutrino is 1.3 x 10 -2 eV, while the mass of the heaviest neutrino is 
3.4 x lO" 2 , 5.8 x 1(T 2 or 9.4 x 1(T 2 eV. 
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1 Introduction 



In the work of Lehmann, Newton, and Wu pQ, the Kobayashi-Maskawa [2j matrix is ex- 
pressed in terms of the masses of the three generations of quarks: 




This is accomplished by introducing a new horizontal symmetry. Some of the earlier 
attempts in this direction are given in [Hill]. 

Recent experiments at Super-Kamiokande [SHZI indicate the presence of neutrino os- 
cillations, which would imply that the neutrinos are not all massless. If it is accepted 
that the neutrinos are not massless, then it is most natural in the Standard Model [Hj to 
introduce three right-handed neutrinos in addition to the three known left-handed ones. 
In this way, there are six quarks and six leptonsQ In this paper, the consequences of a 
universal quark-lepton mixing are studied. In other words, the method of [l a is used to 
express the lepton KM matrix and the neutrino mixing matrix in terms of the masses of 
the three generations of leptons: 




Of course the masses of the three charged leptons are accurately known, leaving as unknown 
parameters the masses of the three neutrinos. Thus there are three parameters to be 
determined instead of seven, the three masses plus the four in the lepton KM matrix. 

It is the purpose of this paper to use the data from solar neutrinos [101 HH H2 IE] and 
atmospheric neutrinos jSJ [7j to determine the three neutrino masses separately, not only 
the differences of their squares. 



2 The mixing matrix 

The quark mixing matrix proposed by Lehmann et al. pQ| relates the mixing to the actual 
quark masses. These mass matrices are each given by four parameters a, b, c, and d, and, 
when applied to the leptons without modification, take the form 

/ d{£) \ 
M{£) = d{£) c{£) b{£) , (2.1) 
\ b{£) a{£)) 

/ id(v) \ 
M{y) = -id{y) c{y) b(v) , (2.2) 
V b(v) a(u)J 

1 As shown in J5J, neutrino-oscillation experiments cannot distinguish between massive Majorana and 
Dirac neutrinos. 
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with 

b 2 (£) = 8c 2 (£), b 2 (u) = 8c». (2.3) 

The diagonalization of these mass matrices is achieved by the orthogonal matrices 
R(£) and R{v), where 

M(£) = R(£)M diag (£)R T (£), (2.4) 

M{u) = diag(-2,l,l) J R(z/)M diag (z/) J R T (z/)diag(^,l,l). (2.5) 

If there is CP violation in the lepton sector, then the imaginary entries in ()2.2j) and (|2.5jl 
are required. If by any chance there is no CP violation in the lepton sector, then both i 
and —i there should be replaced by 1. 

The diagonal mass matrices have the form 

/Ai \ /mi \ 

M diag = A 2 0=0 -m 2 , (2.6) 

\o o A 3 y v° W 

where A 2 < and 

mi < m 3 . (2-7) 

We shall here only be concerned with M(y) and R(y). 

In dealing with quarks pQ, the observed quark masses allow a much stronger form for 
the inequality (|2.7)L namely 

mi < m 2 < TO3. (2-8) 

The lack of direct experimental data on the neutrino masses implies that (|2.7|) can be 
used, but not (|2.8jl . The first task is therefore to determine the allowed region in the space 
(mi,m 2 ,m 3 ), which must be between those permitted by (|2.7|) and (|2.8jl . 

The parameters a, 6, c, and d are related to the masses by the following conditions, 

a + c = Si = m 3 — 1712 + mi, 
8c 2 + d 2 — ac = —5*2 = m 3 m 2 — m 3 mi + rri2mi, 

ad 2 = —S3 = mim^m?,. (2-9) 
The cubic equation for the parameter a is then 

9a 3 - 17.V 2 + (8Sf + S 2 )a - S 3 = 0. (2.10) 

Any real cubic equation can have either one or three real solutions. Where there is one 
real solution, that one is negative, and thus unphysical, as is seen from (|2.9|) . Where there 
are three real solutions, one of them is negative, while two are positive. We shall refer to 
these two positive solutions as Solution 1 (larger a) and Solution 2 (smaller a). 

These considerations can be used to determine the allowed physical region in the 
(7711/7713,777.2/777.3) plane, as shown in Fig. ^ This region is only slightly larger than the 
triangle given by the inequality ()2.8|) . with two additional regions, one where m 2 > m 3 and 
the other a very small one with mi > m 2 . 
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Figure 1: The allowed region for the three neutrino masses mi, m 2 and 7713 is within the 
solid contour. 



3 The three-family MSW mechanism 

The coupled equations satisfied by the three neutrino wave functions are [Hj 





Ml 



M? 2 
M 2 32 



M 2 3 
ML 







fMr)\ 


)} 


1 


Mr) 






W r )J 



(3.1) 



where D(r) = V2GpN e (r), with Gp the Fermi weak-interaction constant and N e (r) the 
solar electron density at a distance r from the center of the sun. Furthermore, we denote 

= 01, ^ = 02, = 03- 

The evolution of the neutrino wave functions is determined by the squared mass matrix, 




b 2 + c 2 + d 2 




'Ml 
M 2 2l 
■ Ml 



Ml 
M 2 2 
Ml 



Ml 
Ml 
ML 



M 2 



= (M 2 



'V 



(3.2) 



The 



the neutrino momentum, p, and the solar electron density. Here, 
eigenvalues of the squared mass matrix (multiplied by r /2p, with r defined below) are 
denoted fix, fi 2 , and /13, and ordered such that 



Mi < \H S ^3- 



< 



(3.3) 
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It is actually a good approximation to take an exponential electron density, N e (r) = 
N e (0) exp(— r/r ). A fit to the solar density as given by jT3] leads to r = 6.983 x 10 4 km. 
For this case of an exponential solar density, the three-component wave equation can be 
solved in terms of generalized hypergeometric functions, 2-^2 |16j . 

We scale and shift the radial variable, u = t/tq + Uo, with uq determined such that 



D(0)r e"° = 1. 
The above equation (|3.1j) may then be written as 



(3.4) 






X-2 




M 


X2 







U>2 


X3 







w 3 



(3.5) 



where uo 2 and u 3 denote the eigenvalues of the 2x2 mass matrix involving ip 2 and ^3. 
With x = e~ u , one finds for ipi a differential equation of the form 



n 



U=i 

the solutions of which, 



dx 



n 



i=2 



dx 



^i = 



(3.6) 



3=1 



(i), 



u 



can be given in terms of generalized hypergeometric functions 2 F 2 as 



1 -i(u 2 - fij), l-i(u 3 -Hj] 
1 — i(jik - ft), 1 — «(^~ ^j) 



le 



(3.7) 



(3. 



with k, £ = 1, 2, 3, k ^ j, £ 7^ j. For the other flavours, ^(u) and ipz{u) are given by 
similar expressions, with shifted parameters. 

There is not much information about these functions 2^2 i n the literature jTTj. For 
the case of two flavours, the products in ()3.6|) go only up to j = 2, and a familiar con- 
fluent hypergeometric function \F± (also denoted Whittaker function or parabolic cylinder 
function) is obtained [T%] . 

In order to impose the boundary conditions that only electron neutrinos are produced in 
the sun, we have to determine these functions at large and negative values of u. The series 
expansion is in principle convergent, but it is not practical for large absolute values of both 
parameters and the argument. Instead, methods have been developed for the evaluation 
by a combination of relating them to 3F1 and asymptotic methods |16j . 
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4 Results 



Since the only unknown parameters for the present theory of universal quark-lepton mixing 
are the three neutrino masses, it remains to determine these masses using experimental 
data on solar and atmospheric neutrinos. 

For the atmospheric neutrinos we take the 8 data points for and the 8 data points 
for v ey as reported recently [Sj. These sixteen data points are treated as separate inputs, 
but we allow two overall normalization constants for the two sets of data. 

For the solar- neutrino data, we use the total rates from the Chlorine experiment [Til] , 
the Gallium experiments fTJ (we average the two results) and the Super-Kamiokande 
experiment [U]. We adopt the neutrino energy spectra and detector efficiencies as given by 
Bahcall et al., and, for the latter detector, we also include the neutral-current cross section 

UJ. 

The determination of the three neutrino masses is to be carried out by a x 2 fit to these 
pieces of data. Since the number of degrees of freedom is 14, a good x 2 fit would give 
support to this idea that the mass mixing is universal for quarks and leptons. For this 
purpose, we have scanned the entire (m 1; m 2 , m 3 ) parameter space. It is the necessity to 
cover this entire space that makes it essential to develop the theory described in Sec. 3 
here. 

It is found that there are three minima. They are: 

(A) Solution 1: (B) Solution 2: (C) Solution 2: 

m 1 = 0.0016 eV m x = 0.0011 eV m x = 0.0011 eV 

m 2 = 0.013 eV m 2 = 0.013 eV m 2 = 0.013 eV 

m 3 = 0.034 eV m 3 = 0.058 eV m 3 = 0.094 eV 

with x 2 = 31 with x 2 = 23 with x 2 = 20 

A few simple conclusions can be drawn from this set of mass values. First, with 14 
degrees of freedom, the three values of x 2 must be considered to be quite good. This 
is evidence in favour of universal quark-lepton mixing. In our mind the difference in 
these three values of x 2 is n °t significant. Secondly, the three values for m 2 are the 
same — a surprise to us. This presumably means that the relative accuracy in the present 
determination of m 2 is better than those of mi and m 3 . This can be understood in the 
following way. The relative accuracy for nil is low because mi itself is quite small, and 
the value of x 2 is n °t very sensitive to such small masses. The average of these three mi 
is given in the abstract, and very roughly we guess the error to be about a factor of two. 
On the other hand, the uncertainly in the value of m 3 is simply due to the three minima 
giving different values. It is a matter of which one of the three minima is physically the 
correct one. 

For completeness, we give the rotation matrix R(u). The Euler angles (8, <p, ip) [20] are: 
(22°, 28°, -8°) for (A), (35°, 21°, -5°) for (B), and (31°, 19°, -3°) for (C). 

To demonstrate the scanning of the entire region of allowed values of neutrino masses, 
we show in Fig. [2 an example of such a scan showing contours of constant x 2 ■ m this 
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figure, the value of m 3 used is 0.058 eV, corresponding to the minimum (B) above. It is 
from such scans that we know there are only three minima. In this figure, some minor 
irregularities along the edges of the allowed region are artifacts of the finite grid spacing. 
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0.04 



0.06 
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Figure 2: Fits to the atmospheric- and solar- neutrino data. Contours are shown at \ 2 — 30, 
35, . . . , 60, 70, 100, 150, . . . , 350. (The outer contour outlines the boundary of the allowed 
region, cf. Fig. [T]) 

It is difficult to compare the present result with the previously given allowed regions 
in the parameter space. The reason is that the allowed regions have typically been given 
on the basis of the mixing of two neutrino species, while there is significant mixing among 
all three neutrinos, or the studies concentrate on either solar or atmospheric neutrinos 
|21[ I22j. The following comparisons are nevertheless of interest. 

(a) It is probably correct to compare the m\ — m\ here with the Sm 2 from atmospheric- 
neutrino data. In this case, the m 2 — m\ for the three minima covers roughly the allowed 
range of 5m 2 , with minimum (B) near the center, and minima (A) and (C) near the edges 
of the allowed region. 

(b) It is probably not too far wrong to identify m 2 — m 2 here with the 5m 2 from 
solar-neutrino data. If so, the values are reasonably close, except with the so-called "low- 
mass-low-probability" (LOW) solution |2*T| . 

(c) It is more difficult to discuss the strength of the coupling. It can nevertheless 
be concluded, in the context of solar neutrinos, that the present solutions are significantly 
closer to the "large-mixing-angle" (LMA) solution than to the "small-mixing-angle" (SMA) 
solution. 
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In connection with the problem of distinguishing between the three minima, more 
accurate data from Super-Kamiokande and related experiments are needed. Also, a better 
understanding of the high-energy hep neutrino flux [21] could make the electron recoil 
energy spectrum useful for this purpose. Furthermore, important information may be 
forthcoming also from the long-baseline experiments and exotic atom experiments 
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